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■ Abstract 

\^ ■ The value of the proton singlet axial constant is estimated. It has been shown 

that the axial anomaly plays a crucial role in this calculation. Obtained result is 
^ , sufficiently suppressed in comparison with the naively expected one. The magnitude 
^ ■ of the strange quark contribution for the proton state is also computed approximately. 

m 



In the last time the problem of calculation of the nucleon singlet axial constant 



r-| ■ has intensively been discussed in the literature. On the classical level, this constant 



determines the value of spin carried by the valence quarks inside the nonrelativistic 

p 1' 

D I nucleon. However, recent experimental data [y lead one to the quantity which is 
much smaller than the naive theoretical value. The authors of ref. (see also 
^ ■ [0]) have ffistly pointed out that it is the axial anomaly which suffers the naive 
^ ■ interpretation at the quantum level. As a result, the apparent "spin crisis" was 



resolved 0]. 

In spite of this, the problem of calculation of the nucleon singlet axial constant 



is far from being solved successfully. Based on the detailed analysis |0|-[|T31) one can 
conclude that both perturbative and nonperturbative contributions are essential for 
this calculation. 

In the present paper, we will try to estimate the value of the proton singlet axial 
constant taking into account characteristic peculiarities of the singlet axial channel. 

Essentially, the problem is as follows: experimental measurements have al- 
lowed one to calculate the integral over the Bjorken variable of the first structure 
function of polarized deep-inelastic lepton-nucleon scattering . In accordance with 
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the Ellis- Jaffe sum rule |Tl]|, this integral is related to the axial constants for the 
corresponding nucleon 

Jo 12 6b 

+C'{l,a,{Q'))^-~gf\ (1) 

Here C^^{l,as{Q^)) and C^{l,as{Q^)) are the nonsinglet and singlet operator co- 
efficient functions normalized to unity at the tree level. The perturbative expansion 
for the nonsinglet coefficient function is known at the next-next-to-leading approxi- 
mation of perturbation theory ||12| whenever it is known at the next-to-leading order 
for the singlet one |jl3|. The axial constants in eq. (1) are normalized as follows: 



< P, 5|43)|p, s >=< P, S\uj^-f,u - h^i,d\P, S >= gfS^, (2a) 



< P, |P, S >=< P, ^|m7^75M + d^,75d - 2s7;.75s|P, S >= gfs^, (26) 



< P, ^|4°)|P, S >=< P, S\u^^^,u + d^^^.d + S7^75S|P, S >= gfS^, 



(2c) 



with Sfj, being the proton spin four-vector . Here the state vectors are normalized in 
the Fock space in accordance with the following condition: 

< P(A;i)|P(A;2) >= {27i)V^6^^\k2 - h)- 

The renormalization group invariant quantity g^^^ in eq. (1) can be written through 
the scale dependent axial constant g^^^ [jl4 



.(0) -r-^^'^dx (0) 

9a =e Jo ^ gy, 
where 7(x) is the singlet axial current anomalous dimension known at the three-loop 



approximation [|I^ and P{x) is the QCD /^-function determined at this order too 
[0. Let us denote the ffist and second addendum in the r.h.s. of eq.(l) as a^"^ and 
, respectively. As a result, Ellis- Jaffe sum rule will take the form 



glix, Q'')dx = a^^(l, a,(Q^)) + a^(l, a.(Q^)). 



(3a) 
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In these notations singlet part of the Elhs-Jaffe sum rule in the leading order appro- 
ximation can be written in the form 

a^ = (l- 0.3333(^i^^))i#. (36) 

TT 9 

Substituting the numerical values for the integral in the l.h.s., for the coefficient 
function C^^{Q'^) and also for the nonsinglet axial constants [|T^, one obtains 
a surprisingly small value for the singlet part ~ 0.02 ± 0.01. It is easy to see 
that the perturbative corrections presented in (3a) do not save the situation and the 
only possibility is to assume (in contrast with the naive expectation) that the proton 
singlet axial constant g\ itself is a small quantity. In what follows we will try to 
estimate the magnitude for ^f^. 

Let us begin with the consideration of the anomaly equation for the flavor singlet 
axial current with three active flavors 

d^Af= Y: 2tm,qj,q + ^GG. (4) 

q=u,d,s 

Brief mention should be made of an important role of the mass terms in this ex- 
pression. At first glance it would seem that these terms are negligibly small due 
to the smallness of the corresponding quark masses. Mass terms are often dropped 
out in the chiral limit. However, it has been demonstrated in [l^ [|l^ |^ that due 
to the fermionic zero modes over the instanton background the matrix elements of 
the operators like q'j^q have singular behavior in the limit of zero quark mass. So, 
only the whole quantity mqq'y^q has a strict sense in the chiral limit; moreover, this 
combination gives a nonzero contribution when the quark masses are set to be equal 
to zero after all calculations are performed. 

On the other hand, one can assume that the mass term is absent in the Lagrangian 
from the very beginning. If so, the derivative of the axial current would be expressed 
through the gluon operator only. However, this is the case when the QCD Lagrangian 
is chiral invariant and consequently the corresponding generating functional could not 
be able to produce a nonzero value for the quark condensate (if the ordinary Feynman 
boundary conditions are assumed). The fulfillment of the axial Ward identities in 
this case becomes doubtful too (for the detailed discussion see ref. ||20|). Another 
way of putting it is that one must calculate quasiaverages instead of simple averages 



2l|] going to the limit rUg ^ after this p2[. Such a procedure allows one to take 



into account consequences of the spontaneous breaking of the chiral symmetry in 
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the pure QCD. It has been demonstrated in |]T9[ that the contributions coming from 
the mass terms in the anomaly equation lead one to the nonrenormalization of the 
^-term in QCD. The mass terms play an important role in our calculation too. 

Now we are in a position to consider the matrix element of the singlet axial 
current over the proton states 

< P{h)\A^^^\P{k,) >= gf{q^)U{k{)^,^,U{h) - tq,gP{q')Uih)^,Uik2), (5) 

where = ki^ — k2^, U{k) is the proton spinor. It is well known that there are no 



massless excitations in the flavor singlet axial channel even in the chiral limit ||23 |. 
Consequently, the phenomenological expression for the pseudoscalar constant gp can 
be written in the form 

^P = E^^, (6) 

where are the mass squares of the corresponding pseudoscalar particles in the 
flavor singlet case. These quantities are not equal to zero even in the chiral limit. 
This is an important point where the singlet and nonsinglet constant calculations 
drastically differ from each other. An analogous expression for the nonsinglet pseudo- 
scalar constant could be presented as a sum of terms containing poles in and the 
terms containing the multiplier of an order of o(m^); as a result, the mixing between 
the axial and the pseudoscalar constants will take place within the corresponding 
sum rules. However, as we will see later, this is not true for the flavor singlet case. 
Reason is that there are no poles in eq. (6) even in the chiral limit. The coefficients 
Bn in this expression are determined by the proton-meson interaction vertices. Using 
now eq. (4), the matrix element (5) takes the form 



< P{k^WAP\P{k,) >= {-zgT{q')2m, + q'gP{q'))U{k^h,U{k: 



(7) 



where rup is the proton mass. Our aim will be to estimate the value for g^ by means 



of the QCD sum rule method ||2^. Following the ideology of this approach consider 
the three-point correlation function 



Tf{k,M) 



Ak\x—i 



< 0\T*7]''{x)A^^\0)f{y)\0 > dxdy 



-iiv < {]\T*7^^{x/2)Af{y)f]^{-x/2)\{] > dxdy, 



(8) 
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where p = {ki + k2)/2, and ri°'{x) is an interpolating current for the proton state |2 

<0|r/"(x)|P(A;i) >= X^U^^ih). 
The phenomenological expression for this correlator can be written in the form 

ih, k,) - (p _ ^2)2 + - (9) 

where k'^ = k\ = k^ = + (f and the kinematical condition (pq)=0 is assumed 
for simplicity. Dots denote here the contributions with the one-pole term and the 
terms without poles in [k"^ — m^. Multiplying expression (9) over iq^ and using (7) 
we obtain 



iq 



^Tf{k,M 



{k? — mp)2 



Let us consider now the operator product expansion (OPE) for our three-point cor- 
relation function. (Hereafter, we will follow the definitions of |^). The bilocal OPE 
for the correlator under investigation looks like 

iq^Tf{k,M)\\p^\^oo = 1 6*^"+^"^ < 0|r*r/"(x/2)9MW(|/)r/'3(-x/2)|0 > dxdy\\,2\^^ 



= Y,CnW I e'^y < 0|T*6„(0)9M(f)(y)|0 > dy+^i?„(p,g) < 0|6„(0)|0 > . (11) 

n n 

The first sum in the r.h.s. of this expansion is determined by the vacuum expectation 
values of some bilocal operators and the second one gives the contributions of the 
local operator vacuum expectation values. The operators O„(0) are generated within 
the OPE for the proton interpolating currents 

I f e'P^T*nx/2)f{-x/2)dx = ECn(p)O„(0) 

It is worthwhile to point out here that the bilocal VEVs play the crucial role in 
this approach, in fact, they give the dominant contribution in the case of the sum 
rules for the proton pseudoscalar constants (for a detailed discussion see ref. [Q). 
However, we will demonstrate that the bilocal VEVs cancel each other in the singlet 
case and only the local part of the OPE determines the value for the proton singlet 
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axial constant. Let us turn to the consideration of the correlators presented in the 
r.h.s. of eq. (11) 

Pniq') = I J e''' < 0|T*O„(0)9M(°)(|/)|0 > dy. 

The dispersion relation for this quantity looks like 

„ / „2 2\ 1 /" IrnPn(s)ds , 

FniQ =—(!) = — 7^ 1" subtractions. 

71 J s + 



Some remarks concerning the peculiarities of this expression are in order here. It is 
well known that the correlator for two singlet axial currents < A^fM^,"^ > contains 
the so-called Kogut-Susskind pole |^ beyond the intermediate physical particle 



contributions. This pole term is produced by the collective excitations over the com- 
plicated vacuum p8|. Within the momentum representation for < A^^M^J') > the 



Kogut-Susskind pole gives an additional term proportional to {q^iqu)/(l'^- Due to sin- 
gularities in g^, the existence of such a contribution makes it difficult to use the sum 
rule formalism. However, in our case, after multiplication of the corresponding corre- 
lator over the Kogut-Susskind contribution becomes the polynomial in momenta 
and, consequently, could be absorbed in the subtractional part of the dispersion re- 
lation. Then, using the Borel transformation this polynomial disappears from the 
sum rules under investigation. Later on, we shall follow this way. 

On the other hand, the phenomenological expression for the imaginary part of 
the correlator under consideration looks like 

a Nf 

ImPn{s) = 7rA(s) < 0|O„|0 >< 0| V 2imqq-f^q + ^^GG\0 > + 

q=u,d,s 

k q=u,d,s 

where A(s) determines the complicated vacuum contribution. The second addendum 
in the r.h.s. of this relation gives the contribution only to the proton pseudoscalar 
constant gp and it is the first term in the r.h.s. which is responsible for g^. However, 
following flSl this last is equal to zero. Indeed, let us integrate over the fermionic 
fields in the path integral representation for < 0|g75g|0 > with the subsequent cal- 
culation of the corresponding determinant in terms of the gluon fields. As an answer 
one can obtain [jTH[ JSg 

as GG 



< 0|g75g|0 >= 



lOTC nig 
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Consequently 

^ 01 ]\fr 

< 0| V 2imqq-i^q + -^-^GG\Q >= 0. 

4:71 



q=u,d,! 

We can conclude that the bilocal operator VEVs do not affect the sum rule for the 
proton singlet axial constant. It is important to point out that this statement is based 
on the assumption of absence of the massless excitations in the singlet channel in the 
chiral limit (relation with the U{1) problem, see ref. [^). Now it is evident that the 
proton singlet axial constant must be small in comparison with the nonsinglet one. 
The point is that the local part of the OPE determining the value for g\ gives the 
contributions of the order only. Let us turn to the consideration of these terms. 
For practical calculations of the local part it is convenient to introduce the following 
quantity 

d 



iXy^{0)2ml^rl5 , C7.75 



where C is some unknown constant. The leading contributions to the sum rule for 
g\ come from the three-loop perturbative diagram and the one-loop diagram with 
the four-fermion condensate. Our main assumption is that the perturbative three- 
loop contribution is suppressed by the loop-factor l/(47r^)^ and could be neglected in 
the case of rough estimation. The only thing we have done is the calculation of the 
four-fermion condensate contributions. In such an approximation, the corresponding 
sum rule takes the form 



Xlg\{Q)2ml iC _ 2 ."^^2 1 < >^ 



where Cp is the Casimir operator of the fundamental representation of SU (3)c group 
[Cf = 4/3) and = -p^. Multiplying this expression over the quantity P^ + mi 
and using then the Borel transformation, we obtain the following sum rule for the 
proton singlet axial constant 

2AJ^° (0)e-^ ^ Cl(^f\ < uu {r- - r), (13) 

where r = M'^/rrip and is the Borel parameter. Substituting the known values 
for the quark condensate and all parameters in eq. (13) |^ we obtain the following 
result 
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^7A(0)U=iGeye-^ ~ 0.025(r2-r). 

The relative stability in these sum rules is reached when r ^ 1.3 — 1.5 and conse- 
quently 

^A(0)U=iGey- 0.02 - 0.03. 

So our crude estimation leads to the value for the proton singlet axial constant which 
is sufficiently small in comparison with unity and the nonsinglet constants. 

Now, using the obtained value for the proton singlet axial constant one can esti- 
mate the s - quark sea contribution to the proton matrix element. Indeed, writing 

one can 



the relation g^^^ — g^^^ = 3 As and substituting the known value for g^^^ [|1^ 
obtain that As|^=iGeV' — —0.16. This estimate is in qualitative agreement with the 
value obtained from the precision fit of the experimental data . 

In conclusion, let us make remark concerning our consideration. The result ob- 
tained here is qualitative in its essence and could be considered as some crude ap- 
proximation only. For a more precise calculation of the value for the proton singlet 
axial constant it is necessary to take into account both the three-loop perturbative 
contribution and the dimension- seven quark-gluon operator contribution within the 
QCD sum rule approach. However, one could not expect that these contributions 
will lead to sufficient increasing in the value for the proton singlet axial constant. 
Missing contributions are of an order of only. We have demonstrated that correct 
treatment of the anomalous contributions leads to sufficient numerical suppression 
for the value of the proton singlet axial constant (in terms of sum rules the can- 
cellation of the bilocal contributions takes place). This statement is based on the 
assumption that there are no massless excitations in the flavor singlet channel even 
in the chiral limit. 

So, we have estimated the value for the proton singlet axial constant. This value 
is sufficiently suppressed in comparison with the naively expected one. The value of 
the s - quark sea contribution for the proton matrix element is also estimated. 

The author is grateful to S.A. Larin, A. A. Pivovarov and A.N. Tavkhelidze for 
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partly supported by the Russian Fund of the Fundamental Research, Grants N 94- 
02-04548a, N 94-02-14428 and "Bazis Bank" graduate student fellowship. 
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